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We theoretically study coherent subharmonic (multi-photon) transitions of a harmonically driven 
spin. We consider two cases: magnetic resonance (MR) with a misaligned, i.e., non-transversal driv¬ 
ing held, and electrically driven spin resonance (EDSR) of an electron conhned in a one-dimensional, 
parabolic quantum dot, subject to Rashba spin-orbit interaction. In the EDSR case, we focus on 
the limit where the orbital level spacing of the quantum dot is the greatest energy scale. Then, 
we apply time-dependent Schrieffer-Wolff perturbation theory to derive a time-dependent effective 
two-level Hamiltonian, allowing to describe both MR and EDSR using the Floquet theory of period¬ 
ically driven two-level systems. In particular, we characterise the fundamental (single-photon) and 
the half-harmonic (two-photon) spin transitions. We demonstrate the appearance of two-photon 
Rabi oscillations, and analytically calculate the fundamental and half-harmonic resonance frequen¬ 
cies and the corresponding Rabi frequencies. For EDSR, we Hnd that both the fundamental and 
the half-harmonic resonance frequency changes upon increasing the strength of the driving elec¬ 
tric field, which is an effect analogous to the Bloch-Siegert shift known from MR. Remarkably, the 
drive-strength dependent correction to the fundamental EDSR resonance frequency has an anoma¬ 
lous, negative sign, in contrast to the corresponding Bloch-Siegert shift in MR which is always 
positive. Our analytical results are supported by numerical simulations, as well as by qualitative 
interpretations for simple limiting cases. 

PACS numbers: 71.70.Ej 73.21.La 76.20.-|-q 


I. INTRODUCTION 

Magnetic resonance (MR) is an established method to 
coherently control the quantum state of spins. A sim¬ 
ple example is a spin-1/2 electron subject to a time- 
dependent magnetic field^^^, described by the Hamilto¬ 
nian 

H =■ cr, (1) 

where the magnetic field B(t) = (0, —i?ac coswt, H) con¬ 
sists of a ‘longitudinal’ static component B and a ‘trans¬ 
verse’ ac component characterised by the drive strength 
Sac a-nd the drive frequency w, and couples to the elec¬ 
tron spin represented by the vector cr = (cr^, cry, cTz) of 
Pauli matrices. 

A typical initial-value problem considered in MR is 
when the initial state of the spin ijj{t = 0) = | f) is 
the ground state of the static Hamiltonian, 
and driving is switched on abruptly at t = 0. In the 
case of weak driving Bac “C B, the rotating wave ap¬ 
proximation (RWA) often provides a satisfactory descrip¬ 
tion of the dynamics. Using this approximation, one 
finds the following simple phenomenology. If the reso¬ 
nance condition huj = ggBB is fulfilled, the drive will 
induce complete Rabi oscillations resulting in a transi¬ 
tion probability = |(j, |'!/(t))p = sin^(Ht/2), where 
Ul = ggBBg,c/(f2h) is called the Rabi frequency. Oth¬ 
erwise, i.e., in the case of a finite detuning d = uj — 
gg,BB/h between the drive frequency and the resonance 
frequency, one finds incomplete Rabi oscillations with a 5- 


dependent frequency: Pi{t) = P™^^ s\v?{y/5“^ -k VlH/2), 
with P'^^^ = OV(02 -b (5^) < 1 . 

Still focusing on the weak-driving regime Bac ^ B, 
one can go beyond the RWA, e.g., by numerical sim¬ 
ulations or analytical techniques such as the Floquet 
perturbation theory^. Then, a richer phenomenology 
is revealed, including (i) subharmonic or ‘multi-photon’ 
resonances^, (ii) drive-strength-dependent Bloch-Siegert 
shifts^ (BSSs) of the resonance frequencies, and (iii) 
Bloch-Siegert oscillations modulating the simple Rabi 
oscillations^. We restrict our attention to (i) and (ii) 
here. 

(i) In the case of a transverse ac field, such as the 
example used in Eq. (I), odd subharmonic resonances 
appear^. Rabi oscillations are obtained not only for 
the fundamental resonance w ~ gg-BB/h, but also when 
u! ~ ggBB/{Nh) with N = 3,5,7,.... In the case of 
a misaligned, non-transversal, ac field, such as B(t) = 
(0, —Bac cos^coswt, B — Bac sinflcos cut) with 0 < 9 < 
TT{2, both even and odd subharmonics are present. The 
Rabi frequency Ul[es at the Al-photon subharmonic res¬ 
onance is weaker than that of the fundamental one: 

oc Bi^/B(^-i). 

(ii) The resonance frequencies uj^es (i-e., the drive fre¬ 
quencies where complete Rabi oscillations are induced) 
increase with increasing drive strength, by an amount 
that depends on TV, and is proportional to Bf^/B. 

In many situations, it is be more convenient to con¬ 
trol spins using an ac electric field rather than an ac 
magnetic field. For example, if an electron spin is elec- 
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trostatically confined in a quantum dot (QD), then an 
ac electric field can be easily created by applying an ac 
voltage component of the confinement gate electrodes. 
Along these lines, electrically driven spin resonance®^^^ 
(EDSR) of individual electron spins was demonstrated in 
a variety of materials^^^^^. As the ac electric field cou¬ 
ples to the orbital degree of freedom of the electron and 
has no direct effect on the spin, a sufficiently strong cou¬ 
pling mechanism between the orbit and spin is required 
for EDSR. Such a coupling can be supplied by spin-orbit 
interaction, hyperfine interaction, or an inhomogeneous 
magnetic field. 

Recent experimental^°’^^“^® and theoretical^^^^^ stud¬ 
ies addressed subharmonic resonances in EDSR. One 
mechanism that leads to subharmonic resonances in 
EDSR is the appearance of higher harmonics Nw of the 
drive frequency ui in the induced orbital dynamics^®’^^. 
In this case, the time-dependent effective magnetic fields 
caused by the orbital dynamics will also have compo¬ 
nents at frequency Nuj^ leading to Rabi oscillations as 
Nuj matches the Zeeman splitting. Higher harmonics in 
the orbital dynamics arise naturally if the confinement 
potential is anharmonic, or if the driving electric field is 
inhomogeneous. Subharmonic EDSR resonances can also 
arise in the presence of harmonic confinement and homo¬ 
geneous ac electric field, if the gradient of the effective 
magnetic fields is inhomogeneous; this is the case, e.g., 
if the the effective magnetic field is spatially localized or 
disordered^^. A third mechanism, able to cause strong 
subharmonics with large is provided by Landau-Zener 
dynamics in the vicinity of level anticrossings^®’^^. 

In this work, we theoretically describe the character¬ 
istics of the half-harmonic resonance in MR with a mis¬ 
aligned ac field, as well as in EDSR. First, we use Flo- 
quet perturbation theory^ to characterise the parame¬ 
ter dependence of the half-harmonic resonance frequency 
and the corresponding Rabi frequency in the case of MR; 
in particular, the BSS is calculated. As for EDSR, we 
study a model®® (see Fig. I), where a single electron is 
parabolically confined in a one-dimensional (ID) quan¬ 
tum dot, and is subject to a dc magnetic field, an ac 
electric field, and spin-orbit interaction of Rashba type, 
the latter three being spatially homogeneous. We show 
that the half-harmonic resonance does arise in this model, 
despite the harmonic confinement and homogeneous ac 
electric field. In the perturbative regime of this model, 
i.e., when the orbital level spacing Hujq dominates over 
other energy scales, we analytically derive the parame¬ 
ter dependence of the half-harmonic resonance frequency 
and the corresponding Rabi frequency. This is achieved 
via a combination of time-dependent Schrieffer-Wolff per¬ 
turbation theory (TDSW), which is used to obtain a 2 x 2 
effective ‘two-level’ or ‘qubit’ Hamiltonian TLq, and Flo- 
quet perturbation theory, applied to describe the qubit 
dynamics governed by Hq. 

This paper is structured as follows. In Sec. H we sum¬ 
marize our main results. Secs. HI and IV are dedicated 
to the detailed discussion of MR and EDSR, respectively. 



FIG. 1. Electrically driven spin resonance in a ID quantum 
dot. (a) An electron occupying the ground state ipoiz) of a 
parabolic confinement potential is excited by an ac electric 
field of amplitude Aac and frequency uj. The electron is sub¬ 
ject to a homogeneous magnetic field B and spin-orbit inter¬ 
action characterised by the direction nso. (b) Left: Orbital 
levels labelled by the oscillator quantum number n, separated 
by the level spacing huiQ. Right: Diagram representing one 
of the many fifth-order virtual processes contributing to the 
half-harmonic resonance. Horizontal lines represent the en¬ 
ergy eigenstates of the harmonic oscillator Hamiltonian. Ar¬ 
rows labelled by E, B and SO correspond to matrix elements 
of the ac electric field, magnetic field and spin-orbit Hamil¬ 
tonians, respectively. Note that the spin-orbit Hamiltonian 
provides both spin-flip and spin-conserving matrix elements 
only if both cos 9 and sinS are nonzero. [Cf. Eq. (36).] 


In them we formulate the problems, derive analytical so¬ 
lutions and compare these with numerical simulations 
where called for. In Sec. V we give a conclusion of our 
findings. 


II. SUMMARY OF THE RESULTS 


In this Section, we summarize the main results that 
are derived in later Sections. Let us start with the case 
of spin-1/2 MR with a misaligned ac field. The param¬ 
eters of the model are B, the energy scale of the static 
magnetic field along the z axis; i?ac, the energy scale of 
the ac magnetic field oriented in the yz plane; and 9, the 
angle enclosed by the ac field and the y axis. (For more 
details, see Sec. HI.) The fundamental (or single-photon) 
resonance frequency Wres, that is, the drive frequency at 
which the Rabi oscillations are complete, deviates from 
the Zeeman splitting: 

^res = R + ^BSS’ (^) 


where the second term on the right hand side is the BSS 
and has the form 


hjj 


( 1 ) _ 
BSS — 


i?ac COS^ 6 

lel 


(3) 
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The Rabi frequency at the fundamental resonance is 


^ cos 9. (4) 

Similarly to the fundamental resonance, the half¬ 
harmonic (two-photon) resonance also acquires a positive 
BSS: 


where 


B 


f^iel = ^ + 


BSS 


(5) 

( 6 ) 


The Rabi frequency at the half-harmonic resonance is 


hnill 


Rfc sin 2d 
Ib 


(7) 


Note that the resonance frequencies above are expressed 
up to second order in the small energy scale Sac *C B. 
For a detailed discussion of these results, see Sec. HID. 

In the case of EDSR in a ID parabolic QD, the param¬ 
eters characterizing the model are as follows. The orbital 
level spacing fcjg is the dominant energy scale; the static 
magnetic field, oriented along the z axis, is character¬ 
ized by the Zeeman-splitting energy scale B; spin-orbit 
interaction is described by the energy scale a and the 
unit vector ngo = (0, cosd, sind) which points along the 
spin-orbit field; and the energy scale E^c describing the 
strength of the driving ac electric field. We find the fol¬ 
lowing results for the fundamental resonance frequency: 

=B + /iwW -p -b (8) 

where the correction consists of a g-factor renormaliza¬ 
tion term. 


, fi-) 2Bd^cos^d/ a^(l-b sin^ d)\ 

=-- —) 

a term describing the non-linear Zeeman effect. 




2B^a^ cos^ d 


(9) 


( 10 ) 


and a correction that is second order in the drive strength 
Bac, hence analogous to the BSS: 


( 1 ) Ba^El^cos^e 

- ^4^4 


( 11 ) 


Note that the sign of this BSS is negative, in contrast to 
the positive sign in the case of MR [Eqs. (3) and (6)]; an 
interpretation of this anomalous sign is given in Sec. V, 
paragraph (4). The Rabi frequency at the fundamental 
resonance is 


hOS}} = 2 


BEscOl cos d 


1-b 


B 2 - 252 \ 
h?ujl j ■ 


( 12 ) 


The half-harmonic resonance frequency is shifted with re¬ 
spect to the half of the fundamental resonance frequency: 


= \{^ + + ^niz) + ^BSS> (13) 

where the drive-strength-dependent BSS is expressed as 

(14) 


(2) 2Ba'^E'^ cos^ d 


ftWBSS — 




The Rabi frequency at the half-harmonic resonance is 


(2) _ Bd^B,^,sin(2d) 

(ftwo)4 ■ ^ ^ 

Note that the EDSR resonance and Rabi frequencies 
above are expressed up to fifth order in the small energy 
scales a, B,Bac ^o- A. detailed discussion of these 
results, and their comparison with numerical solutions of 
the time-dependent Schrddinger equation, is included in 
Sec. IV D. 


III. MAGNETIC RESONANCE WITH A 
MISALIGNED AC FIELD 

In this Section, using Eloquet perturbation theory, 
we derive and discuss the properties of the fundamen¬ 
tal (single-photon) and half-harmonic (two-photon) res¬ 
onances in MR, for the spin-1/2 case. In particular, the 
results (2) — (7) are derived. 


A. Problem formulation 

We consider MR spin dynamics driven by a misaligned 
ac field. The Hamiltonian reads 

B{t) =■ cr , (16) 

where the magnetic field has the form 

B(t) = I — BacCOSdcOSWt . (17) 

\ B — Bac sin d cos ojt j 

Here we introduced B = g^sB and Bac = gpisBac. 
Henceforth the parameters with tilde, e.g., B, have en¬ 
ergy dimension. The Hamiltonian in Eq. (16) has 
four parameters: the strength of the static field B, the 
strength of the driving field Bac, the frequency of the 
driving field w, and the misalignment angle d. Note that 
d = 0 corresponds to a transverse ac field, and d = 7r/2 
corresponds to a longitudinal ac field. We consider the 
case of weak driving, Bac ^ B. 

In particular, we want to solve the initial-value prob¬ 
lem described in Sec. I: the initial state is the ground 
state I t) of the Hamiltonian without driving i.e., '0(^ = 
0) = I t); driving is switched on abruptly at t = 0, 
and we are interested in the time evolution ^/>(t) of this 
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state. We calculate the transition probability describing 
the time-dependent occupation of the excited state | J,) 
at the fundamental and half-harmonic resonances, and 
from those we deduce the parameter dependencies of the 
resonance frequencies and the Rabi frequencies. 

In the rest of this Section, we use Floquet perturba¬ 
tion theory^’^®^^® to derive the results and to provide 
qualitative interpretations in simple limiting cases, such 
as the limits of transversal and longitudinal ac fields. 
Even though similar treatments can be found in the 
literature^’®®, we present a detailed discussion of the MR 
problem for the following reason. The MR problem is rel¬ 
atively simple as compared to the EDSR problem, which 
can be appreciated, e.g., by comparing the driven two- 
level Hamiltonians of Eqs. (16) and (41), respectively. 
Moreover, as we will show, the Floquet method and the 
qualitative interpretations we describe here for the MR 
problem can be carried over to the EDSR problem, once 
the 2x2 effective qubit Hamiltonian in Eq. (41) is ob¬ 
tained for the latter. This allows us to provide a rather 
compact description of the EDSR in the forthcoming Sec¬ 
tions, by referencing this Section wherever possible. 

B. Floquet method 

The Floquet method allows one to find the solution 
of an initial-value problem of a periodically driven quan¬ 
tum system, described by the time-periodic Hamiltonian 
H{t) = H{t + T). The period of the driving is de¬ 
noted by T, and the corresponding (angular) frequency 
by a; = 27r/T. The key ingredient of the method is the 
quantum-mechanical Floquet theorem®®, which guaran¬ 
tees that the Schrodinger equation ih^{t) = 'H{t)'^{t) 
of a d-level system has d solutions ^'^(t) (fc = 1,..., d) 
that are themselves periodic with period T, apart from a 
phase factor. Therefore, these special solutions have the 
form 

d oo 

|«'fe(t)) ^ , (18) 

m— — oo 

where \tjji) is an arbitrary basis of the Hilbert space. Note 
that the result of the double sum is a periodic function 
of t with period T. In Eq. (18), the quantity Ek and 
the coefficients Ck,im are a priori unknown; the former is 
called quasi-energy. Once these special solutions |'I'fc(f)) 
are found, they provide the propagator 

d 

c/(t,0) = ^|vkfe(t))(vl/,(0)|, (19) 

which in turn provides the solution of any initial-value 
problem via 

4'(t) = [/(t,0)4'(0). (20) 

The special solutions 'kfc(t) are found by using Eq. (18) 
as an Ansatz, substituting it to the Schrodinger equation. 


evaluating the scalar product of the equation with (^;/1, 
multiplying the equation by e”™ and integrating the 
equation in time between t = 0 and t = T. This proce¬ 
dure yields the following eigenvalue equation for Ek- 

d oo 

E E ,lmCk,lm — EkCk,l' m' 1 (21) 

1 — 1 m— — co 

where 

OO 

'Fl'm' ,lm — TTlhwdl/iSfjiiYTi E 

n— — oo 

is the Floquet matrix or Floquet Hamiltonian, and we 
introduced the Fourier components of the Hamilto¬ 
nian via 

OO 

nit) = (23) 

n— — oo 

We call two eigenvalue-eigenvector pairs of F equiva¬ 
lent, if the two time-dependent solutions they generate 
via Eq. (18) are the same. Importantly, even though the 
number of eigenvalue-eigenvector pairs of T is infinite, 
they form only d equivalence classes. 

In summary, we have transformed the time-dependent 
Schrodinger equation of the periodically driven d x d 
Hamiltonian 'H(t) into the time-independent Schrodinger 
equation (21) of the infinite-dimensional Floquet Hamil¬ 
tonian F. To construct the special solutions (18), and 
thereby the solution of any initial-value problem via Eqs. 
(19) and (20), the quasi-energies Ek and the correspond¬ 
ing eigenvectors Ck should be found by solving the eigen¬ 
value problem of the Floquet Hamiltonian F. 

C. Perturbative description of the transition 
probability 

After reviewing the Floquet method in general, we now 
apply this to the MR problem defined in Eq. (16). Here, 
we have a two-level system, therefore d = 2, and we use 
|a) = I t) = Ifpi) and 1,5) = | ),) = \tfj 2 ) to denote these 
levels. The Fourier components of the Hamiltonian read 

^ (24a) 

= ii?ac(cosdtTy -|-sind(Tz), (24b) 

and the other Fourier components are zero. 

First, consider the case when the drive frequency is 
close to the fundamental resonance, fua k, B. Then, the 
diagonal elements of the Floquet Hamiltonian F {Floquet 
levels) form pairs: 

Fam,am — TflhlO ^ —1- 

(25) 

The distance between different pairs is approximately 
fijjj K, B, which is much larger than the energy scale Hac 
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characterising the off-diagonal elements of J-. Therefore, 
the tools of quantum-mechanical perturbation theory can 
be used to provide an approximate solution of the eigen¬ 
value problem of the Floquet Hamiltonian. 

The structure of the Floquet Hamiltonian T is visu¬ 
alised for the case fujj = B in the level diagram shown 
in Fig. 2. Horizontal lines represent the diagonal ma¬ 
trix elements !Fim,im of the Floquet Hamiltonian, their 
vertical positions correspond to their value, their colour 
(black, red) represents their spin index I G (a, /?) and 
their horizontal position stands for their Floquet index 
m = ..., — 1,0,1, 2,.... The vertical spacing of the Flo¬ 
quet levels is hu! = B. The blue arrows indicate the 
nonzero off-diagonal matrix elements of which are of 
the order of Hac and hence small compared to the level 
spacing. 

In the case huj = B shown in Fig. 2, the Floquet lev¬ 
els form degenerate pairs. The pair formed by 
and J'a,o,a,o is highlighted in Fig. 2 by the blue box. The 
subspace of this pair is weakly coupled to the other Flo¬ 
quet levels, hence this coupling can be treated perturba- 
tively using (time-independent) Schrieffer-Wolff pertur¬ 
bation theory^®, which is also known as quasi-degenerate 
perturbation theoryThis perturbative treatment is 
also applicable if there is a finite, but small detuning 


8 = io — B/h <C B/h from the resonance condition. 
The small dimensionless parameter characterising the 
strength of the perturbation is e = B^^c/B. 

In this case, the Floquet Hamiltonian reads 


Fundamental resonance 



FIG. 2. Magnetic resonance in a misaligned ac field: struc¬ 
ture of the Floquet Hamiltonian at the fundamental reso¬ 
nance. Panels show cases when the ac field is perpendicular 
to the static field (a), is parallel to the static field (b), has 
finite perpendicular and parallel components (c). Horizon¬ 
tal lines (blue arrows) correspond to diagonal (off-diagonal) 
matrix elements of the Floquet Hamiltonian B. The vertical 
position of each horizontal line corresponds to the value of 
the diagonal matrix element. 


J 




a_i 

P-i 

ao 

/3o 

ai 

Pi 




i 





a_i 


. . . —\B — hu! 

0 

ii?acsin0 

-|i?ac COs6» 

0 

0 

/3-1 


0 

^B — hw 

\Bs.cCOs9 

-iHacSin6> 

0 

0 

ao 


... ii?acsin6> 

-|.BacCOs6< 

-iB 

0 

iHac sin0 

-|i?acCOs6» ... 

Po 


... |HacCOs6» 

-ii?acSin0 

0 

\B 

fHac cos 9 

-iHacSin0 ... 

ai 


0 

0 

ii?acsin0 

-|i?ac C0S6» 

—\b + hio 

0 

/3i 


0 

0 

|HacCOs6< 

-ii3acSin6> 

0 

hw 


(26) 


1. Fundamental resonance within RWA 


Using first-order perturbation theory, the two non¬ 
equivalent eigenvalues and eigenvectors of T can be found 
approximately. For this we introduce Bq and so that 
J- = J-Q + J- 1 . is the diagonal component of T at 

uj = B/h., i.e. at (5 = 0. 


First-order perturbation theory in J-i amounts to di¬ 
agonalizing the 2x2 block highlighted in purple in Eq. 
(26) and Fig. 2. For future reference, we recast this 2x2 
block to the form 


where eo 


^_ cq -f A iX 

—iX eo — A ’ 

— ^{B + hS), A = —h6/2, and A 


(27) 

ii?acCOs6>. 
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The matrix T has eigenvalues 

S± = eo ± \/A2 + A2 , (28) 

and corresponding eigenvectors 

£± = A±[^(A±\/A2 + A2),1], (29) 

where is a normalization constant. Note that in¬ 
stead of using the numerical index k G (1,2) labelling 
the solutions (18), in Eqs (28) and after we use the val¬ 
ues k G (-I-, — ). 

The results (28) and (29) imply that the two non¬ 
equivalent approximate eigenvalue-eigenvector pairs of 
are {E±,c±), where 

( c±,i if (/,m) = (/3,-l), 
c±,im = < c ±,2 if (/,m) = (a,0), , (30) 

I 0 otherwise. 

and c± 1 and c± 2 are the components of c± in Eq. (29). 
This result allows us to construct the transition probabil¬ 
ity P/ 3 <_a(t) = |(/3|'I'(t))p from the initial spin (ground) 
state I t) = |a) to the excited state | |) = |/1) via Eqs. 
(18), (19) and (20). A straightforward calculation yields 

= ^^^sin2 . (31) 

According to Eq. (31), the spin makes complete Rabi 
oscillations if A = 0, that is, S = uj — B jh = 0. Hence the 
single-photon resonance frequency is Wres = B/h. The 
Rabi frequency upon resonant driving is fires = 2\/h = 
cos 6] thus only the transverse component of the ac 
field contributes to the Rabi frequency at the fundamen¬ 
tal resonance. In fact, the result (31) is equivalent to the 
one obtained by neglecting the longitudinal ac field and 
performing RWA. 


2. Fundamental resonance beyond the RWA: Bloch-Siegert 
shift of the resonance frequency 


Let us discuss the corrections to Wres and fires beyond 
the RWA. To this end, we incorporate in the analysis 
the effect of those matrix elements of Pi that connect 
the two highlighted Floquet levels [see Eq. (26) and Fig. 
(2)] to the complementary subspace. This is done via a 
(time-independent) Schrieffer-Wolff transformation that 
is second order in Pi. The resulting effective 2x2 Floquet 
Hamiltonian P has the form given in Eq. (27), with 


A = - 



cos^ Q 
32H 


A = 


Hac COS Q 

4 


(32a) 

(32b) 


Recall that the eigenvalues and eigenvectors of P are 
given by Eqs. (28) and (29). From these, we conclude 


that the two non-equivalent approximate eigenvalue- 
eigenvector pairs of P are {E±,c±), where 

f c±,i + o(e2) if (/,m) = (/3,-1), 
c±,im = < c ±,2 + o(e^) if (/,m) = (a, 0), (33) 

[ o(e) otherwise. 

We neglect the perturbative corrections ~ o(e),o(e^) in 
the eigenvectors c±, and this implies that the approxi¬ 
mate transition probability is given by Eq. (31). Equa¬ 
tion (31) predicts that complete Rabi oscillations are in¬ 
duced when A = 0; solving Eq. (32a) for uj (recall that 
S = UJ — B/h) provides the resonance frequency shown 
in Eq. (2). The second term of Eq. (2) corresponds to 
the Bloch-Siegert shift of the resonance frequency: as the 
drive strength Hac is increased, the resonance frequency 
shifts upwards. This feature is further discussed in Sec. 
HID. Finally, the Rabi frequency at the fundamental res¬ 
onance, which is given by = 2A, is expressed using 
Eq. (32b) in Eq. (4); the result is the same as in the 
RWA. 


3. Half-harmonic resonance 

Let us now consider the spin dynamics at half¬ 
harmonic resonance, when Huj ~ H/2. The level di¬ 
agram visualising the Floquet Hamiltonian in the case 
huj = B12 is shown in Fig. 3. Again, we can iden¬ 
tify degenerate pairs of Floquet levels, e.g., the pair 
[Pp-I p-i^Pai ai) highlighted with the blue box in 
Fig.’ 3.’ ’ 


Half-harmonic resonance 



FIG. 3. Magnetic resonance in a misaligned ac field: struc¬ 
ture of the Floquet Hamiltonian at the half-harmonic reso¬ 
nance. Panels show cases when the ac held is perpendicular 
to the static held (a), is parallel to the static held (b), has 
hnite perpendicular and parallel components (c). Horizon¬ 
tal lines (blue arrows) correspond to diagonal (off-diagonal) 
matrix elements of the Floquet Hamiltonian P. The vertical 
position of each horizontal line corresponds to the value of 
the diagonal matrix element. 

Note that in this case, there is no direct matrix ele¬ 
ment (blue arrow) connecting these two Floquet levels. 
This implies that by repeating the first-order perturba¬ 
tion theory (equivalent to RWA) done in Sec. HIC 1, we 
would conclude that two-photon Rabi oscillations do not 
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happen. However, this result is not correct; two-photon 
Rabi oscillations can happen. To see that, we perform 
a second-order Schrieffer-Wolff transformation on J", as 
done in Sec. Ill C 2. Furthermore we use the appropriate 
notation S = w — ^. The obtained effective 2x2 Floquet 
Hamiltonian T has the form given in Eq. (27), with 


A J- r ^ac COS^ 9 

A = ns - 


A = 


6B 
sin 29 
85 ■ 


(34a) 

(34b) 


Following the approach used in Sec. HIC2, and solv¬ 
ing A = 0 we find that the half-harmonic resonance fre- 
quency Wres is given by Eq. (5). Furthermore, using 
hQrel = 2A and Eq. (34b), the Rabi frequency at the 
half-harmonic resonance is obtained as shown in Eq. (7). 


D. Discussion 

Let us now discuss the main features of the results (2), 
(4), (5), and (7). 

Consider first the fundamental resonance frequency 
Wres expressed in Eq. (2). The second term in Eq. (2) 
implies that tUres has a positive drive-strength-dependent 
correction oc B^^/B with respect to the nominal Zeeman 
splitting B. This correction is known as the BSS, which 
can be regarded as a special case of the ac Stark shift^^. 

Note that the parameter A and hence the Rabi fre¬ 
quency Hres sets the frequency broadening of the funda¬ 
mental transition, as indicated by the prefactor A^/(A^ -I- 
A^) on the right hand side of Eq. (31). According to Eq. 
(4), this power broadening of the fundamental resonance 
is greater by a factor of BjBg^c than the BSS. 

Equation (2) also shows that the BSS is finite in the 
limit of purely transversal drive (9 = 0), and vanishes 
in the limit of purely longitudinal drive (9 = 7r/2). The 
respective Floquet level diagrams in Fig. 2a and b pro¬ 
vide a straightforward interpretation: the BSS can be 
regarded as a consequence of coupling-induced repulsion 
between the Floquet levels. In Fig. 2a (9 = 0), the 
Floquet level is connected by a blue arrow 

(off-diagonal matrix elements of T) to the lower-lying 
Floquet level J^a,- 2 ,a,- 2 - The consequence of this cou¬ 
pling in second-order perturbation theory is level repul¬ 
sion; i.e., the lower-lying Floquet level pushes 
upwards. Similarly, Tao,aO is pushed downwards by its 
coupling to the higher-lying Floquet level These 

second-order level shifts appear in Eq. (32a) as the last 
term, and give rise to a finite BSS. In contrast, each of the 
highlighted Floquet levels in Fig. 2b {9 = 7r/2) is con¬ 
nected to one higher-lying and one lower-lying Floquet 
level, and the corresponding downward and upward level 
repulsions cancel each other, giving rise to a vanishing 
BSS in this case. 


Consider now the half-harmonic resonance. Equation 
(7) provides the corresponding Rabi frequency, and it in¬ 
dicates the existence of Rabi oscillations unless 0 = 0 or 
9 = 7r/2. I.e., Rabi oscillations appear at half-harmonic 
excitation only if the transversal and longitudinal compo¬ 
nents of the driving field are both nonzero. The Floquet 
level diagrams shown in Fig. 3 provide a visual interpre¬ 
tation of this feature: Rabi oscillations arise if the blue 
arrows (off-diagonal matrix elements of B) draw at least 
one path between the two Floquet levels highlighted by 
the purple box, via virtual intermediate Floquet levels 
outside the box. In the special cases 0 = 0 and 0 = 7r/2 
depicted in Figs. 3a and b, respectively, no such paths 
exist. However, there exist infinitely many such paths for 
0 < 0 < 7r/2 (Fig. 3c), due to the coexistence of spin- 
conserving and spin-flip off-diagonal matrix elements. In 
particular, in our second-order Schrieffer-Wolff transfor¬ 
mation leading to the result (7), the two two-step paths 
via J^a,Q,a,Q and are incorporated. 

In the case of the half-harmonic resonance, the re¬ 
lation between the power broadening and the BSS is 
qualitatively different from the case of the fundamental 
resonance. For the half-harmonic resonance, the power 
broadening is given by Eq. (7), whereas the BSS is given 
by the second term of Eq. (5), i.e., the two quantities 
are of the same order, both being ~ E^^jB. Hence we 
expect that for the half-harmonic resonance, the BSS is 
relatively easily resolvable experimentally, at least if the 
dissipative frequency scales are smaller than the power 
broadening. 

Equation (5) also shows that the BSS is finite in 
the limit of purely transversal excitation (0 = 0), and 
vanishes in the limit of purely longitudinal excitation 
(0 = 7r/2). An interpretation completely analogous to 
the case of the fundamental resonance can be given based 
on the Floquet level diagrams in Fig. 3a and b. 


IV. ELECTRICALLY DRIVEN SPIN 
RESONANCE 

A. The model 


From now on, we describe EDSR mediated by spin- 
orbit interaction in a ID parabolic quantum dot. The 
setup is shown in Fig. 1. The Hamiltonian 


'H=%o+%e+'Hb+'Hso (35) 


includes the harmonic-oscillator Hamiltonian ('Ho) con¬ 
sisting of the kinetic energy of the electron and the 
parabolic confinement potential, the ac electric potential 
arising from the driving electric field {'He), the static 
Zeeman effect caused by a homogeneous magnetic field 
(Hb), and the spin-orbit term (Hso)- The explicit forms 




B. Effective qubit Hamiltonian 


of these terms, respectively, are as follows: 

■Ho = ^ + = huo , (36a) 

He = e^-Eac sin(a;t) = ifac sin(a;t)(a^ + a ), (36b) 
Hb = = -^Ba^ (36c) 

Hso = ctPzi^so • cr = za(a^ — a)nso • cr (36d) 

Here, a and are the ladder operators of the harmonic 
oscillator Hamiltonian, and ngo = (0,cosd,sind) is the 
direction of the effective magnetic field arising from spin- 
orbit coupling. Furthermore, we defined 

B = g*PBB 

mhuiQ 
2 

These quantitites have the dimension of energy. 

Note that we use the same notation 6 for two different 
quantities: 9 appears in Eq. (17) as the ac field misalign¬ 
ment angle in MR, and it also appears in this Section 
and in Fig. 1, as the angle characterising the direction of 
the spin-orbit term. We use the same notation for these 
quantities as they play very similar roles in the spin dy¬ 
namics. 

It is natural to represent the Hamiltonian terms (36) 
in the product basis of the orbital and spin degrees of 
freedom, {\na) | n = 0,1, 2,...; cr =f,4,}, where n is the 
harmonic-oscillator orbital quantum number and a is the 
spin quantum number with quantization along z. 

We will refer to the two lowest-energy eigenenstates 
of our static Hamiltonian Hq + Hb + Hso as the qubit 
basis states. The qubit basis state with the lower (higher) 
energy will be denoted by |G) {\E)). 

The electron is initialized in state |G) at t = 0. Our 
aim is to describe the time evolution of the state upon 
driving. In particular, we are interested in the time- 
dependent occupation probability PE{t) of state \E). It 
is expected that at resonant driving hui se B, the dy¬ 
namics resembles Rabi oscillations. Subharmonic (multi¬ 
photon or A^-photon) resonances at /Ej ~ B/N (where 
N = 1,2,...) are also expected. In this work we focus 
on the fundamental (single-photon, = I) and half¬ 
harmonic (two-photon, N = 2) resonances. 

We aim at an analytical, perturbative description of 
spin transitions induced by the ac electric field. In partic¬ 
ular, we calculate the resonance frequency and the Rabi 
frequency at resonant driving. We consider the parame¬ 
ter range where the energy scale Hujo of the confinement 
potential dominates the other four energy scales, the lat¬ 
ter ones being assumed to be comparable in magnitude: 

fcj ~ d ~ Eac ~ H ^ hwQ. (40) 



(37) 

(38) 

(39) 


This hierarchy of energy scales will allow for a perturba¬ 
tive description of the dynamics, with the small parame¬ 
ter e ^ ^ jr- ^ ^ 1- 

tjjQ fiUjQ hUjQ hUjQ 


In the EDSR problem defined in Sec. IV A, the hier¬ 
archy of the energy scales is given by Eq. (40). Because 
of this hierarchy, an effective time-dependent two-level 
Hamiltonian [see Eq. (41) below] can be derived for the 
qubit dynamics, using TDSW perturbation theory, which 
we outline in Appendix A. This qubit Hamiltonian can 
then be used to express the resonance frequencies Wres 
and Wres, and the corresponding Rabi frequencies at these 
resonances, flrll and Hres, corresponding to the funda¬ 
mental and half-harmonic resonances, respectively [see 
Eqs. (8), (8), (13), and (15) below]. 

We use the orbital-spin product basis {\na) \ n = 
0,1, 2,...; a =t,4.}, as the starting point of TDSW, and 
take the two-dimensional subspace of [O-i-) and [Oj,) as the 
relevant subspace in TDSW. We carry out a fifth-order 
TDSW (in the small parameter e), which is expected to 
describe both the fundamental and the half-harmonic res¬ 
onances. The TDSW procedure yields the effective qubit 
Hamiltonian 


iig ^ ) + 77^3) + H[^^ + H ^^, (41) 


where the six terms, representing terms from different 
orders in the perturbation, are listed below in Eq. (42). 
Note that the terms H^q ^, H^q '^, and H^q^ are proportional 
to the 2x2 unit matrix (Tq, therefore they do not influ¬ 
ence the dynamics, and hence we disregard them in the 
forthcoming calculations; nevertheless we include them 
here for completeness: 


B^q'’ 

77 ( 3 ) 


77W 

Hf'> 


—ao 

(42a) 

B 

(42b) 

+ Ei, sin^iujt) 

hujQ 

(42c) 

BEacCt cos 9 , , 

(42d) 

a cos 9 ~ 

2 (Aachw cos(wt) -1- Ba sin 9)a2 


ta 2 cos^ 9 — Ea.cbio sm.9 cos{wt)) 

n c<Jq 

0-3 

{Bacos9)^ 

(42e) 


(42f) 
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In Eq. (42f), we used 

BaEp,c cos 9 


,(5) Eacdw^COSd 

/.n = ^-^cosM) 


(43a) 


i?Qf^ sin20, ~T T ~o 9/ . 

(-B - sin (wi)) (43b) 


= 


Eacdw^ sin0 

&UQ 

cos^ 6 


cos{ujt) 


—(B -d +B^^sin (wt)). (43c) 


~ f3) 

Note that the upper index in, e.g., 'Hq refers to the order 
of perturbation theory in which the term appears. 

Out of the six terms in Eq. (41), H^q^ and 'H^q'^ are sim¬ 
ply the projected parts of 'Hq and 'Hi = He + Hb+Hso^ 
respectively. Hq contains a static and a time-dependent 
second-order energy shift, due to the spin-orbit inter- 
action and the ac electric field, respectively. Hq has 
five terms. The first, second and fifth terms are spin- 
and time-dependent, hence these all contribute to the 
qubit dynamics. The third and fourth terms are static; 
they describe the spin-orbit-induced g-tensor renormal¬ 
ization. The fourth-order term H^^ of the qubit Hamil¬ 
tonian, being diagonal, does not influence spin dynam¬ 
ics. The static parts of the fifth-order term 'Hq^^ de¬ 
scribe higher-order g-tensor renormalisation (those pro¬ 
portional to a^B), or nonlinear Zeeman splitting (those 
proportional to 

Already at this point, there are reasons to expect 
that in this EDSR model, a half-harmonic resonance 
occurs, and that the half-harmonic resonance frequency 
is driving-strength dependent: (i) The third-order effec- 

tive Hamiltonian Hq ' incorporates both longitudinal and 
transverse ac components, in analogy with the case of 
the misaligned-field MR discussed in Sec. HI. (ii) The 
fifth-order effective Hamiltonian Hq ^ incorporates terms 
proportional to uit = ~ cos2ujt). The 

longitudinal static part oc E^^az can be interpreted as a 
drive-strength-dependent effective g-tensor renormalisa¬ 
tion, which contributes to the BSS, whereas the dynami¬ 
cal part oc cos 2ujt ay is expected to drive Rabi oscil¬ 
lations at half-harmonic excitation, i.e., when 2huj B. 

We note that the effective Hamiltonian in Eq. (42) 
fulfills the expectation that no spin transition occurs if 
the external B-field and the spin-orbit field are aligned, 
ie, when 9 = tx 12. 


C. Floquet perturbation theory for EDSR 

We apply Floquet perturbation theory, outlined in Sec. 
HI C, to describe the fundamental and half-harmonic res¬ 


onances. In particular, we derive the parameter depen¬ 
dence of the corresponding resonance frequencies Wris 
and Wres, as well as the Rabi frequencies flrll and Ores, at 
these two resonances, up to terms of the order of ^ Be'^. 
There are two significant differences in the derivation of 
the EDSR results with respect to that of the MR results; 
we outline these differences in the following. 

(1) The MR Hamiltonian (16) has a driving term that 
is proportional to sinwt. In contrast, the effective qubit 
Hamiltonian (41) we obtained for EDSR has cos ujt terms 
as well as second-harmonic terms proportional to cos 2ujt. 
In practice, the latter fact implies that the Floquet ma¬ 
trix will contain off-diagonal matrix elements that con¬ 
nect Floquet levels with next-nearest-neighbor Floquet 
quantum numbers. 

(2) In the EDSR case, we repeat the same second-order 
time-indepedent Schrieffer-Wolff transformation on the 
Floquet Hamiltonian E that we applied in Secs. HI C 2 
and HIC3. The Floquet Hamiltonian itself contains 
terms of the order of B, Be^ and Be^, since it is 
constructed from the effective qubit Hamiltonian that is 
itself the result of a finite-order perturbative calculation. 
When we separate the Floquet Hamiltonian to diagonal 
(Bo) and off-diagonal (Bi) components, and apply time- 
independent Schrieffer-Wolff transformation up to second 
order in Bi, the resulting 2x2 effective Floquet Hamilto¬ 
nian will involve higher-order terms, up to Be®. As our 
original Hamiltonian was accurate only up to the ^ Be^ 
terms, we drop the terms that are of higher order than 
~ Be^ from the effective Floquet Hamiltonian. 


D. Analytical vs. numerical solution 

The results we obtain from Floquet perturbation the¬ 
ory are shown in Sec. H as Eqs. (8)-(I5). In the rest 
of this subsection we discuss these results and compare 
them to numerical results. 

The terms describing the fundamental resonance fre¬ 
quency in Eq. (8) are interpreted as nominal Zeeman 
splitting, g-tensor renormalisation, nonlinear Zeeman ef¬ 
fect, and BSS, respectively. We call the last term a BSS 
as it is a power-dependent correction to the resonance 
frequency, that is second order in the drive amplitude, 
hence analogous to the BSS in MR. Remarkably, the 
BSS in Eq. (8) is a negative correction, whereas the 
BSS in MR is always positive. The last term of the half¬ 
harmonic resonance frequency [Eq. (13)] is also inter¬ 
preted as a BSS. Further similarities with the MR case: 
(i) For the fundamental resonance, the BSS is smaller 
(~ Be^) than the power broadening, the latter being 
given by ftHres ~ Be^. (ii) For the half-harmonic reso¬ 
nance, the BSS, being ^ Be'^, is comparable to the power 
broadening, the latter being given by /iHres ~ Be^. (iii) 
The BSS for both the fundamental and the half-harmonic 
resonance is proportional to cos^ 0, i.e., it vanishes in 
the limit of purely longitudinal excitation, and finite for 
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FIG. 4. (color online) Electrically driven Rabi spin dynam¬ 
ics at half-harmonic resonance. The excited-state occupation 
probability Pe is shown as a function of the drive frequency 
u! and time t\ the numerical data reveals the chevron pattern 
characteristic of magnetic resonance. Parameters: 9 = 7r/4, 
a/B = Es^jB = 1, huiojB = 5. The analytical results for 
the half-harmonic resonance frequency a;)es and the Rabi fre- 
quency flies are also displayed. For the above parameter val¬ 
ues, the latter one is related to the time period of the oscilla¬ 
tion via 27r/nres = 27r X 625 h/B . 


purely transversal excitation. These features can be ex¬ 
plained by the argument provided in Sec. HID for the 
case of MR, applied to the effective qubit Hamiltonian 
(41). 

Regarding the results (8) and (13) for the resonance 
frequencies, we note that their ratio is exactly two in the 

limit of vanishing driving power, i.e., liin^j ) = 

2 . 

We have checked that the result (8) for the fundamen¬ 
tal Rabi frequency Hres matches the corresponding result 
of Ref. 5; see Appendix B for details. 

The analytical results are tested against numerically 
exact solutions of the time-dependent Schrddinger equa¬ 
tion defined by the Hamiltonian 'H in Eq. (35). The nu¬ 
merical results were obtained using the truncated Hilbert 
space spanned by the 8 lowest-energy eigenstates of 
'Ho +'Hb, corresponding to the 4 lowest-lying levels of 
the harmonic oscillator. We have checked that there was 
no visible change in the numerical results upon extending 
the Hilbert space with further, higher-lying orbitals. 

In Fig. 4, we plot the numerically computed time 
evolution of the occupation probability of the excited 
state \E), for a finite range of the driving frequency in 
the vicinity of the ‘nominal’ half-harmonic resonance fre¬ 
quency huj/B = 0.5 (see caption for parameter values). 
The analytical result (13) predicts complete Rabi oscil¬ 
lations at w = Wres = 0.4809i?/?i, and the Rabi fre¬ 
quency at this resonance is predicted by Eq. (15) to 



Eac/B 

FIG. 5. (color online) Bloch-Siegert shift and power broad¬ 
ening of the half-harmonic resonance. The maximal excited- 
state occupation probability is shown as a function of 

the amplitude E^c of the driving ac electric field and drive 
frequency ui. Parameters: 9 = 7r/4 a/B = 1, and hoJo/B = 8. 
The red line indicates the analytical result for the resonance 
frequency as the function of electric field based on Eq. (13). 

be Ores ~ ■ These predictions are in line with the 

numerical data shown in Fig. 4. For a finite detuning 
from the resonance frequency, the Rabi oscillations be¬ 
come faster and reduced (i.e., they do not reach Pe = 1), 
leading to the characteristic chevron pattern^’^^ known 
from MR. The results of Fig. 4 therefore reveal simple 
Rabi dynamics at the half-harmonic resonance. 

The density plot of Fig. 5 is a visual demonstration of 
the BSS, i.e., of that the resonance frequency increases 
with increasing drive strength. The figure shows the 
maximum of the excited-state probability PE{t) 

within a time span exceeding the Rabi period at the half¬ 
harmonic resonance, as a function of the amplitude Fac of 
the ac electric field and the drive frequency tu. (See cap¬ 
tion for parameters.) Therefore, vertical cuts of the den¬ 
sity plot correspond to resonance curves. The solid line 
represents the analytical result (13) for the half-harmonic 
resonance frequency. The agreement between the ana¬ 
lytical curve and the Ri 1 ridge of the numerical 

simulation reassures the validity and correspondence of 
the two approaches. Importantly, in Fig. 5, the BSS is 
comparable in magnitude to the power broadening, which 
makes the BSS relatively easily resolvable in experiments 
realizing the model we use. 

A further question is how the BSS depends on the an¬ 
gle 9 characterizing the direction of the spin-orbit inter¬ 
action. This dependence is exemplified by Fig. 6, which 
shows as a function of 6 and the drive frequency. 

The latter is measured from half of the calculated fun¬ 
damental resonance frequency Wres, see Eq. (8). The 
solid line, showing good agreement with the centre of the 
bright p™^’' r; 1 region of the underlying density plot, 
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0 7c/8 7i/ 4 37t/8 7 i/2 

0 

FIG. 6. (color online) Anisotropy of the half-harmonic res¬ 
onance. The maximal excited-state occupation probability 
pmax jg shown as a function of the angle 6 characterizing 
the spin-orbit interaction and the drive frequency lo measured 
from uJtII/2, see Eq. (8). The solid line corresponds to the 
analytically obtained half-harmonic resonance frequency, i.e., 

f'Ol (1) ~ ~ 

it shows w^es — see Eq. (13). Parameters: E/B = 1, 
a/B = 1, and huJojE = 5. 


shows the analytical result for the half-harmonic reso¬ 
nance frequency Wres [Eq. (13)]. 

V. DISCUSSION AND CONCLUSIONS 

(1) We provide a numerical example to estimate orders 
of magnitudes of the EDSR resonance shifts and Rabi 
frequencies. Let us take R = 1.7 T with the electronic g 
factor 2, yielding R « 0.1 meV. We set d = 0.1 meV and 
Eac =0.1 meV, and the orbital level spacing is chosen 
to be Hloq = 1 meV. Then the order of magnitude of the 
Rabi frequency at the fundamental resonance becomes 
fires ^ Be^/H Ri 1.5 X 10®i Corresponding to a spin- 
flip time of « 4.3 ns. For the half-harmonic resonance, 
fires ~ Be^/h m 1.5 X lO^i, corresponding to a spin- 
flip time of « 430 ns. For both resonances, the BSS is 
comparable to the value of fires estimated above. 

(2) The results presented in this work describe a per¬ 
turbative regime where spin-orbit interaction is assumed 
to be ‘weak’, in the sense that the spin-orbit energy 
scale in the QD is dominated by the QD level spacing, 
a hujQ. In nanowire QD host materials such as InAs^® 
and InSb^^, spin-orbit interaction is known to be ‘strong’ 
in the sense that it creates a strong g-factor renormali¬ 
sation, already in the bulk materials. A question arising 
from these facts is: are typical InAs and InSb nanowire 
QDs within the range of validity of our perturbative the¬ 
ory? One way to answer this question is via a compar¬ 
ison of the dependence of the fundamental EDSR res¬ 


onance frequency obtained from the perturbative the¬ 
ory and from experiments. The experiments^^’^® have 
found that the fundamental resonance frequency shows 
a similar angular dependence as the perturbative result 
Eq. (8), i.e., for a magnetic field with a fixed magni¬ 
tude, the resonance frequency is maximal if the magnetic 
field is aligned along a certain direction and minimal if 
it is aligned perpendicular to that direction. To be spe¬ 
cific, we take the data given in the first row of Table 
I. of Ref. 25, which indicates that the ratio of the mini¬ 
mal and maximal resonance frequencies in the considered 
case were ~ 0.84. Using the first two terms in Eq. (8), 
we can identify that ratio with 1 — 2a^/h?u}Q, yielding 
ajhioo ~ 0.28 for this particular InAs device. A similar 
analysis of the experimental data in Fig. 3c of Ref. 17 
results in an estimate ajhioo ~ 0.37 for the measured 
InSb device. These estimates suggest that the InAs and 
InSb QDs are on the border between ‘weak’ and ‘strong’ 
spin-orbit interaction. 


(3) To our knowledge, three experiments have re¬ 
ported subharmonic EDSR resonances in semiconductor 
nanowire QDs, where our model based on the Rashba- 
type spin-orbit interaction could be appropriate to de¬ 
scribe the spin dynamics. The strong subharmonic reso¬ 
nances reported in Stehlik et al.^® are described by a the¬ 
ory developed for strongly driven double quantum dots^^. 
Faint half-harmonic resonances are visible in the data 
of Refs. 25 (see Fig. 2b therein) and 17 (see Fig. 2b 
therein). A quantitative experimental analysis exploring 
the parameter dependencies of the corresponding reso¬ 
nance and Rabi frequencies would allow for a comparison 
with our predictions. 

(4) One of our conclusions was that the BSS of the 
fundamental EDSR resonance frequency has an anoma¬ 
lous, negative sign, see Eq. (8). Here, we provide a 
simple physical picture explaining this result, using the 
unitary transformation applied in Ref. 42. For simplic¬ 
ity, we focus on the case when the spin-orbit field is 
perpendicular to the magnetic field, i.e., 9 = 0. Then, 
the unitary transformation S of Ref. 42 (not to be con¬ 
fused with the generator of the Schrieffer-Wolff transfor¬ 
mation in Appendix A) applied on our static Hamilto¬ 
nian Ho +Hb + Hso eliminates the spin-orbit term and 
transforms the homogeneous magnetic field Hb to an in¬ 
homogeneous, spiral-like magnetic field, H'^ = SHbS^ oc 
B [az cos{z/^) — ax sin(z/^)], where ^ oc 1/5 is the spin- 
orbit length [see Eq. (2) of Ref. 42]. The driving elec¬ 
tric field, incorporated in our model as He , induces a 
spatial oscillation z{t) = — Asinwt of the electron’s cen¬ 
tre of mass with an amplitude A oc Eac- Inserting this 
time-dependent z{t) to the above expression for and 
expanding the terms up to second order in A/^, we find 

(X B (72 ^1 — ^sin^wt^ -I-(Ta, j sinwt . That is, 
the time-averaged z component of the time-dependent 
magnetic field in H'g{t) acquires a correction propor¬ 


tional to —B^ oc —BE^^cP'. Notice that this correction 
is negative and has the same parameter dependence as 
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the BSS in, the last term of, Eq. (8). 

(5) To our knowledge, BSS has not yet been experimen¬ 
tally or theoretically analyzed in the context of EDSR. 
However, we wish to point out that certain numerical 
results in Ref. 10, related to EDSR in a double quan¬ 
tum dot, are reminiscent of the BSS. Figures 4a and 4b 
in Ref. 10 show spin Rabi oscillations for different drive 
strengths. Therein, the drive strength is characterised 
by the dimensionless quantity /. In Figs. 4a and 4b of 
Ref. 10 it is shown that the complete Rabi oscillations 
at the fundamental resonance become incomplete upon 
increasing the drive strength from / = 0.02 to / = 0.15, 
while the driving frequency is maintained. This phe¬ 
nomenology is reminiscent of the effect of BSS: when the 
drive strength is increased, BSS provides a shift of the 
resonance frequency, hence a fixed drive frequency be¬ 
comes off-resonant, and the Rabi-oscillation amplitude 
decreases. It is therefore tempting to interpret these re¬ 
sults as consequences of BSS. However, the phenomenol¬ 
ogy of BSS would imply that (i) upon further increase of 
the drive strength, e.g., at / = 0.35, the amplitude of the 
Rabi oscillation further decreases, and (ii) the Rabi os¬ 
cillation speeds up gradually as / is increased from 0.02 
to 0.15 and to 0.35. The results shown in Figs. 4a and 
4b of Ref. 10 disagree with these expectations, hence we 
conclude that the BSS phenomenology is insufficient to 
describe the numerical results of Ref. 10. Importantly, 
Ref. 10 considers parameter settings where the undriven 
system consists of four, approximately equidistant levels 
(see their Fig. 2a), which is a key difference with respect 
to the effectively two-level setup considered in our present 
work, and can be responsible for the phenomenology de¬ 
viating from that of the BSS. 

(6) Even though EDSR experiments can be performed 
on single QDs^^’^^, many current experiments use the 
Pauli blockade setup for initialization and readout^. The 
latter setup consists of a double QD which is occupied by 
two electrons (in the simplest case) during EDSR. The in¬ 
herent anharmonicity of the double QD confinement po¬ 
tential, as well as the presence of the Coulomb interaction 
between the two (or more) electrons, can provide alter¬ 
native nonlinear EDSR mechanisms^’^®’^^’^®’^®’^^, which 
compete with those presented in our work focusing on 
harmonic confinement and single-electron dynamics. For 
example, an apparently well understood^'*’ case when the 
two-electron and double-QD features dominate the sub¬ 
harmonic EDSR resonances is the experiment of Ref. 26. 

In conclusion, we have studied the characteristics of 
EDSR in a ID QD model with parabolic confinement, 
homogeneous Rashba spin-orbit interaction and homoge¬ 
neous driving electric field. We demonstrated the exis¬ 
tence of subharmonic (multi-photon) resonances in this 
model, and analysed the half-harmonic (two-photon) res¬ 
onance in detail. We have analytically described the pa¬ 
rameter dependence of the fundamental resonance fre¬ 
quency and the half-harmonic resonance frequency, and 
demonstrated that these resonance frequencies increase 
with increasing drive strength. This effect is analogous 


to the BSS in MR. 

Our results describe a perturbative regime, where the 
orbital level spacing of the QD dominates the energy 
scales of the external magnetic field, spin-orbit interac¬ 
tion, and electrical drive. Therefore our results have di¬ 
rect experimental relevance for QDs with weak spin-orbit 
interaction. They can also serve as benchmarks for nu¬ 
merical studies departing from the perturbative regime. 
The model used here contains only minimal ingredients 
necessary to describe EDSR, suggesting that the subhar¬ 
monic resonances and the BSS discussed here are generic 
features of electrically driven spin dynamics. 

Note added: Upon completing this work, we became 
aware of a related experiments*^ revealing half-harmonic 
EDSR in a single QD, mediated by an inhomogeneous 
magnetic field. 
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Appendix A: Time-dependent SchriefFer-Wolff 
perturbation theory 

Here we introduce the time-dependent SchriefFer-Wolff 
perturbation theory (TDSW), the method we use to de¬ 
rive the effective 2x2 time-dependent Hamiltonian (41) 
governing the dynamics of the qubit. 

Let us first recall the basic idea of standard time- 
independent SchriefFer-Wolff (SW) perturbation the- 
ory.39do consider a Hamiltonian TL = "Ho + Td' where 
"Ho is diagonal and 'H' is the perturbation. Furthermore, 
the basis states of H can be divided into relevant A and 
irrelevant B subspaces that have well separated energy 
scales. A and B are weakly interacting, i.e. the matrix 
elements connecting them are small compared to the en¬ 
ergy separation of the two subspaces. Ideally, we can 
introduce a unitary transformation e~^ that brings T-L 
into a block diagonal form T-L = e~^'He^ where the rele¬ 
vant and irrelevant subspaces are separated as illustrated 
in Fig. 7. However, in most of the cases we don’t know 
the explicit form of the transformation e~^ so we have 
to construct it bit-by-bit until the elements connecting 
the two subsets vanish up to the desired order of per¬ 
turbation. This is usually done by expanding e~^ in a 
series and constructing the terms of different orders suc¬ 
cessively. 
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FIG. 7. Schematic form of the matrices we encounter during 
the SW perturbation theory. The label e refers to blocks 
with matrix elements that are much smaller than the energy 
difference between the relevant and irrelevant subspaces. In 
the EDSR problem, the energy scale of the those blocks is 
B ehujo. 

A great advantage of SW with respect to conventional 
perturbation theory is that here we don’t need to distin¬ 
guish between the degenerate and non-degenerate cases. 

Now we introduce the time-dependent SW pertur¬ 
bation theory as a natural extension of the time- 
independent case. Similar approaches have been applied 
for particular problems in Refs. 44-46; here, we provide 
a general description of the method, which we utilised in 
the main text for deriving the effective qubit Hamiltonian 
(41) of the EDSR problem. 

Consider the time-dependent Hamiltonian H(t) = 
T-Lo 'H'{t), where the perturbation is divided into a 
block-diagonal and block-off-diagonal part 'H'(t) = 'Hi -I- 
H 2 as shown in Fig 7. 

In our problem (see Sec. IVA), Hi = H_b and 
H 2 (t) = V-Eit) + Hso- Note that there Hi happens to 


be a time-independent perturbation, but the treatment 
outlined here is readily applicable to a time-dependent 
block-diagonal perturbation as well. 

Similarly to the SW we successively build the unitary 
transformation U{t) = that separates the sub¬ 

spaces A and B, but here the matrix S{t) is now time- 
dependent. Note that any unitary transformation can 
be written in this form, and the matrix S{t) should be 
anti-Hermitian to ensure the unitary character of U(t). 
The matrix S{t) is chosen to be block-off-diagonal (see 
Fig 7). Note also that because of the weakness of the 
inter-subspace coupling, the unitary tranformation U{t) 
is close to unity, and hence S{t) is small and can be 
expressed as a power series with respect the perturbing 
terms. 

The transformation of time-dependent Schrodinger 
equation —j/i^'0(t)-|-H(t)^(t) = 0 with the above t/(t) is 
canonical, i.e., it preserves the form of the time evolution 
equation. The transformed wave function and Hamilto¬ 
nian read as: 

m = (Al) 


■H{t) = -b (A2) 

From now on, we might suppress the time argument and 
denote time derivatives such as as ip. 

Starting from Eq. (A2), we utilize the power series of 
the exponential function. The second term in Eq. (A2) 
is the heart of the time-dependent SW transformation; 
in the time-independent case this term vanishes as S is 
time-independent. The expansion of the first term in 
Eq. (A2) is known from SW formalism, therefore we do 
not discuss it here. The explicit form of the second term, 
after expanding the exponential function, has the follow¬ 
ing form: 


de-^ 

dt 




+ .. ) = {-S+lsS+^SS-^SS^-^SSS-^S^S +...) X 




00 


1 


3!' 


i=o 


U + 1)! 


[^,5] 


U) 


(A3) 


The transformed Hamiltonian then equals to 

^ \n,sr'- in y 

I ^' 


00 ^ 00 -| 

« = i; 7T - “V TTXWI®’ ■ (*“> 


1=0 ' 




with = [[H,S']^”\S'] and [H,S']^°^ = H. 

Note that the second term in (A4) is new with respect 


to time-independent SW, and it is a consequence of the 
time dependence of the Hamiltonian and therefore that 
of the matrix S. Considering a time-independent Hamil¬ 
tonian the second term vanishes and we are left with the 
well-known SW transformation. 

We now exploit the block-off-diagonal property of S 
in order to separate the block-off-diagonal and block- 




























14 


diagonal parts of the transformed Hamiltonian: 

1 r„, „l(2i + l) 


OO 

"^off-diag — ^ ^ 
i-0 

OO 


(2j + l)! 


[Hq + HuSY 


3=0 

OO 

-ihY, 


1 


3=0 


(2j + 1)! 


[5,5] 


(2i) 


(AS) 


"Wdiag = ^ ['Ho + Hi, 5] 


(2l) 


1=0 

OO 






1 


1=0 


(2j + 2)! 


[5,5] 


( 21 + 1 ) 


Then, 5 is determined by solving 

Hoff-diag — 0 ■ 


(A6) 


(A7) 


The effective (now block-diagonal) Hamiltonian becomes 
H = Hdiag- Note that H as well as the term ‘effective 
Hamiltonian’ is also used to describe the block of H cor¬ 
responding to the relevant subpsace. 

So far no approximation has been made; now we make 
use of the smallness of the perturbation. Following the 
approach of time-independent SW perturbation theory, 
we aim at solving Eq. (A7) via expanding 5 as a power 
series in the perturbation, 

5 = 5i + 52 + 53 + ... , (A8) 

where Sj represents an operator of jth order in the per¬ 
turbation. Recall that in TDSW, 5 is time dependent, 
and its time derivative appears in its defining equation 
(A7) as well as in the effective Hamiltonian (A6). There¬ 
fore, to separate the terms of different order in perturbing 
parameter in Eq. (A7), it is necessary to make an a pri¬ 
ori assumption on the frequency scale characterizing the 
magnitude of Sj. As the drive frequency is w, expectedly 
the frequency characterizing the time evolution of all Sj- 
s will be ~ w, hence we assume Sj ~ LuSj. In the EDSR 
problem defined in Sec. IV A, the relevant subspace is 
the subspace of the ground-state orbital spanned by [Of) 
and [Oj,). Furthermore, the energy scales of the drive fre¬ 
quency, drive strength, Zeeman splitting and spin-orbit 
coupling are much lower than the splitting between the 
oscillator levels ~ Wg, and all of them are treated as per¬ 
turbation. This implies that Sj is of the order of (j -I- 1) 
in perturbation. 

Obviously, after solving Eq. (A7) with this assump¬ 
tion, we need to check if the obtained Sj functions are 
consistent with our assumption above. 


From the order-by-order expansion of Eq. (A7), we ob¬ 
tain the following hierarchy of simple algebraic equations 


for the Sj matrices: 



[Ho, 5i] = 

-H 2 , 


(A9a) 

[Ho, 52 ] = 

— [Hi, 5i] -I- 

ihSi , 

(A9b) 

[Ho, 53 ] = 

-[Hi,52]- 

^[n2,Sif'> + tHS2,{^^c) 

[Ho, 54 ] = 

-[Hi,53]- 

- [[H 2 ,5i], 52 ] 



-g [[H2,52] 

, 5i] -1- ihSs , 

(A9d) 


Once the first equation (A9a) is solved for 5i(t), the 
solution can be inserted to (A9a) which then forms an 
algebraic equation for 52 (t), etc. Note that since we work 
in the eigenbasis of Ho, the above procedure simplifies 
to subsequently solving single-variable linear equations, 
which is a trivial analytical task, well suited for symbolic 
computation. 

After obtaining the Sj matrices and inserting them 
into Eq. (A6), we have an order-by-order expansion H = 


Hdiag = H^”\ where 

= Ho (AlOa) 

■H(^)=Hi (AlOb) 

H^') = [H 2 , 5i] + i [Ho, (AlOc) 

H(3) = [n 2 , 52] + i [Hi, 5i]('^ + i [[Ho, 5i], 52] 

+i[[Ho,52],5i]-z;ii[5i,5i] (AlOd) 


With the use of Eqs. (A9) we can further simplify 


Eqs. (AlO): 




h(°) 

= Ho 



(Alla) 


= Hi 



(Allb) 

^(2) 

= ^ [H2,5i] 



(Allc) 

^(3) 

= ^ [H2,52] 

1 

" 4! 


(Alld) 

■^(4) 

= ^[H2,53]- 

[H2,5i](3) 

(Alle) 

■^(5) 

= ^[H2,54]- 

1 

" 4! 

[[[H2,5i],5i],52] 

(Allf) 


[[[H2, 5i] , 52], 5i] - i [[[H2,52], 5i], 5i] 
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Finally, we need to check the consistency of our as¬ 
sumption for the time evolution of Sj with the actual so¬ 
lution we obtained for Sj using that assumption. From 
(A9a), Si inherits harmonic time-dependence from 'H 2 
with frequency uj. This implies that the time derivative 
is ~ ujSi, as assumed. From Eq. (A9b), the matrix 
S '2 might contain frequency components at w, as well as 
at zero frequency and 2a; (if 'Hi is time-dependent with 
frequency uj); nevertheless, the S 2 ^ UJS 2 relation still 
holds, etc. 

In conclusion, TDSW allows for obtaining an effective 
time-dependent Hamiltonian for the relevant subspace. 
The procedure is to evaluate the transformation matrices 
Sj up to the desired order via solving Eq. (A9), and 
substituting the resulting Sj matrices into Eq. (All). 

Appendix B: Rabi frequency of the fundamental 
resonance: relation to the results of Ref. 5 


y — z plane, and the E-held is along the x axis. Further¬ 
more, we project the Hamiltonian on the j/-ground-state 
orbital of the harmonic oscillator, yielding 

p 2 1 2 

Hgbl = ^ + + aPxCTy -f • cr 

+ eEs_cX shoot (Bl) 

For simplicity, we focus on the special case B = (0,0, B) 
from now on. Then, the Hamiltonian in Eq. (Bl) is 
equivalent to our Hamiltonian H at 0 = 0. 

To deduce the Rabi frequency calculated by GBL for 
the special case above, we start from their Eqs. (13) and 
(14), where they provide the time-dependent part of the 
effective qubit Hamiltonian as 

where 


EDSR in a QD in a two-dimensional electron gas due to 
Rashba and Dresselhaus spin-orbit interactions has been 
described by Golovach, Borhani and Loss (GBL) in Ref. 
5. Therein, the Rabi frequency of the fundamental reso¬ 
nance as a function of system parameters (magnetic held 
strength, magnetic held direction, spin-orbit interaction 
strengths and ac electric held amplitude and direction) 
has been calculated. Even though the dimensionality and 
the spin-orbit Hamiltonian in the model of GBL differs 
from our model, the calculated Rabi frequencies can be 
compared after a special case of the model of GBL has 
been reduced to one dimension. Here we show that after 
this dimension reduction our result for the fundamental 
Rabi frequency equals that of GBL. 

In the model of GBL, the 2DEG lies in the x-y plane. 
We consider the special case when the conhnement po¬ 
tential is parabolic and has a cylindrical symmetry, the 
Dresselhaus coupling vanishes, /? = 0, the B-held is in the 


h(t) = 2/rBB X (B3) 

A straightforward calculation shows that 

aeE^cg*y,BB . 

-h{t) = -—^- sm{u}t)e^ (B4) 

2 iMjJq 

= 2^^^ sm{ujt)ea; (B5) 

Note that this effective ac magnetic field is perpendicular 
to the static magnetic field, which is applied in the z 
direction. The Rabi frequency due to this ac magnetic 
field at the fundamental resonance frequency reads 


n 


( 1 ) 

res,GBL 


h^UJQ 


(B6) 


which is identical to our result in Eq. (8), if the latter 
is evaluated at 0 = 0 and terms above third order are 
dropped. 


^ F. Bloch and A. Siegert, Phys. Rev. 57, 522 (1940). 

^ J. H. Shirley, Phys. Rev. 138, B979 (1965). 

® F. H. L. Koppens, C. Buizert, K. J. TielrooiJ, I. T. Vink, 

K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 

L. M. K. Vandersypen, Nature 442, 766 (2006). 

^ M. Veldhorst, J. C. C. Hwang, C. H. Yang, A. W. Leenstra, 
B. de Ronde, J. P. Dehollain, J. T. Muhonen, F. E. Hudson, 
K. M. Itoh, A. Morello, and A. S. Dzurak, ArXiv e-prints 
(2014), arXiv: 1407.1950 [cond-mat.mes-hall]. 

® V. N. Golovach, M. Borhani, and D. Loss, 
Phys. Rev. B 74, 165319 (2006). 

® C. Flindt, A. S. Sprensen, and K. Flensberg, 
Phys. Rev. Lett. 97, 240501 (2006). 

^ Y. Tokura, W. G. van der Wiel, T. Obata, and S. Tarucha, 
Phys. Rev. Lett. 96, 047202 (2006). 

® E. 1. Rashba, Phys. Rev. B 78, 195302 (2008). 

^ J. D. Walls, Phys. Rev. B 76, 195307 (2007). 

D. V. Khomitsky, L. V. Gulyaev, and E. Y. Sherman, 


Phys. Rev. B 85, 125312 (2012). 

“ R. Li, J. Q. You, G. P. Sun, and F. Nori, 
Phys. Rev. Lett. Ill, 086805 (2013). 

Y. Kato, R. C. Myers, D. G. Driscoll, A. G. Gossard, 

J. Levy, and D. D. Awschalom, Science 299, 1201 (2003). 

K. G. Nowack, F. H. L. Koppens, Y. V. Nazarov, and 

L. M. K. Vandersypen, Science 318, 1430 (2007). 

M. Pioro-Ladriere, T. Obata, Y. Tokura, Y.-S. Shin, 
T. Kubo, K. Yoshida, T. Taniyama, and S. Tarucha, 
Nat. Phys. 4, 776 (2008). 

E. A. Laird, C. Barthel, E. 1. Rashba, G. M. 
Marcus, M. P. Hanson, and A. G. Gossard, 
Phys. Rev. Lett. 99, 246601 (2007). 

S. Nadj-Perge, S. M. Frolov, E. P. A. M. Bakkers, and 

L. P. Kouwenhoven, Nature 468, 1084 (2010). 

S. Nadj-Perge, V. S. Pribiag, J. W. G. van den 
Berg, K. Zuo, S. R. Plissard, E. P. A. M. 
Bakkers, S. M. Frolov, and L. P. Kouwenhoven, 






16 


Phys. Rev. Lett. 108, 166801 (2012). 

F. Pei, E. A. Laird, G. A. Steele, and L. P. Kouwenhoven, 
Nat. Nanotech. 7, 630 (2012). 

V. S. Pribiag, S. Nadj-Perge, S. M. Frolov, J. W. G. 
van den Berg, 1. van Weperen, S. R. Plissard, 
E. P. A. M. Bakkers, and L. P. Kouwenhoven, 

Nat Nano 8, 170 (2013). 

E. A. Laird, F. Pei, and L. P. Kouwenhoven, 

Nat. Nanotech. 8, 565 (2013). 

E. Kawakami, P. Scarlino, D. R. Ward, F. R. Braakman, 

D. E. Savage, M. G. Lagally, M. Friesen, S. N. Cop¬ 
persmith, M. A. Eriksson, and L. M. K. Vandersypen, 
Nat. Nanotech 9, 666 (2014). 

P. V. Klimov, A. L. Falk, B. B. Buckley, and D. D. 
Awschalom, Phys. Rev. Lett. 112, 087601 (2014). 

F. Forster, M. Miihlbacher, D. Schuh, W. Wegscheider, 
and S. Ludwig, Phys. Rev. B 91, 195417 (2015). 

E. A. Laird, C. Barthel, E. 1. Rashba, C. M. Marcus, M. P. 
Hanson, and A. C. Gossard, Semiconductor Science and 
Technology 24, 064004 (2009). 

M. D. Schroer, K. D. Petersson, M. Jung, and J. R. Petta, 
Phys. Rev. Lett. 107, 176811 (2011). 

J. Stehlik, M. D. Schroer, M. Z. Maialle, M. H. Degani, 
and J. R. Petta, Phys. Rev. Lett. 112, 227601 (2014). 

A. De, C. E. Pryor, and M. E. Flatte, 
Phys. Rev. Lett. 102, 017603 (2009). 

J. Pingenot, C. E. Pryor, and M. E. Flatte, 
Phys. Rev. B 84, 195403 (2011). 

E. 1. Rashba, Phys. Rev. B 84, 241305 (2011). 

M. P. Nowak, B. Szafran, and F. M. Peeters, 

Phys. Rev. B 86, 125428 (2012). 

E. N. Osika, B. Szafran, and M. P. Nowak, 


Phys. Rev. B 88, 165302 (2013). 

G. Szechenyi and A. Palyi, 

Phys. Rev. B 89, 115409 (2014). 

E. N. Osika, A. Mrerica, and B. Szafran, 

Phys. Rev. B 90, 125302 (2014). 

J. Danon and M. S. Rudner, 

Phys. Rev. Lett. 113, 247002 (2014). 

M. Trif, V. N. Golovach, and D. Loss, 

Phys. Rev. B 77, 045434 (2008). 

T. Dittrich, Quantum Transport and Dissipation (Wiley- 
WCH, 1998). 

P. K. Aravind and J. O. Hirschfelder, 

The Journal of Physical Chemistry 88, 4788 (1984). 

I. Gromov and A. Schweiger, J. Magn. Res. 146, 110 

( 2000 ). 

J. R. Schrieffer and P. A. Wolff, 

Phys. Rev. 149, 491 (1966). 

R. Winkler, Spin-orbit Coupling Effects in Two-Dimensional Electron 
Springer Tracts in Modern Physics, Vol. 191 (Springer, 

2003). 

C. Wei, A. S. M. Windsor, and N. B. Manson, 

J. Phys. B: At. Mol. Opt. Phys. 30, 4877 (1997). 

L. S. Levitov and E. 1. Rashba, 

Phys. Rev. B 67, 115324 (2003). 

P. Scarlino, E. Kawakami, D. R. Ward, D. E. Savage, 

M. G. Lagally, M. Friesen, S. N. Coppersmith, M. A. Eriks¬ 
son, and L. M. K. Vandersypen, ArXiv: 1504.06436 (un¬ 
published). 

F. Schwabl, Advanced Quantum Mechanics, Chapter 9. 

(Springer, 2008). 

A. Kaminski, Y. V. Nazarov, and L. 1. Glazman, 

Phys. Rev. B 62, 8154 (2000). 

Y. Goldin and Y. Avishai, Phys. Rev. B 61, 16750 (2000). 


32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 


